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Abstract 

Effect of external periodic force on an oscillatory order in a reaction diffusion 
system (Gierer Meinhardt model) has been investigated. The 2:1 resonance situa- 
tion is found susceptible for the generation of a band of phase instabilities. These 
phase instabilities, captured on multiple time scales, produces a mismatch between 
the oscillation frequency of reacting species. 

PACS number(s): 87.10.+e, 47.70.Fw 

The effect of parametric periodic forcing on oscillatory reaction diffusion systems are 
being studied with renewed interest to see frequency entrainment and resulting multi- 
phase, steady as well as traveling, orders separated by phase fronts [1, 2, 3, 4]. Existence 
of multiphase oscillations are theoretically accounted for by showing the stability of phase 
separated oscillatory orders in complex Ginzburg Landau equation or in some reaction 
diffusion models. The mechanism that can possibly cause a slow drift in overall phases 
of oscillation under periodic forcing, and thus produce stable phase-separated regions, is 
an important subject for investigation. In view of that, we are going to investigate the 
effect of periodic forcing in time, on an oscillatory system, on multiple time scales. In al- 
most all reactions diffusion systems, one of the reacting species is dependent on the other 
for its production and thus does not need be externally supplied. This situation causes 
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a constant phase difference to appear between the homogeneous oscillatory reactants of 
the system. Here we focus on the possible slow variation of that phase difference as a 
consequence of varied response to the applied force by the reacting species. The reaction 
diffusion system we have worked on, is the Gierer-Meinhardt model. The results obtained 
are the generation of a band of phase instabilities caused by parametric resonance at 2:1 
sub-harmonic response to the applied frequency. The instabilities developed at slow time 
scales and produce a mismatch in frequencies of oscillation of the activator and inhibitor 
species. The result is interesting because in such a situation when activator and inhibitors 
are oscillating with different frequencies, many different things can follow. An incommen- 
surency in the activator and inhibitor frequencies of oscillation can possibly generate weak 
phase chaos. In other case, if one of the two species tries to lock itself in phase with the 
other, a slow drift in the overall phase of oscillation can presumably be an outcome. In 
order to get rid of such frequency mismatch at various orders, a finer adjustment of the 
parameter is needed. Such an adjustment also causes a nontrivial shift in phase boundary 
obtained from the linear stability analysis of the system. 

The Gierer-Meinhardt model that we have taken up is [5, 6] 



D is the diffusivity of the activator A and is always less than unity to satisfy the Turing 
condition. The a is the basic production rate of activator where /x can be interpreted 
as the production constant and at the same time removal rate of the inhibitor B. The 
linear stability analysis of this model shows that the homogeneous stationary basic fixed 
point B — — [1 -\- aY which becomes unstable [5] to a time independent spatially 
inhomogeneous state when 



dA 
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and to a time periodic spatially homogeneous state if 



l-a 



1 + a 



(3) 



For the time periodic state the oscillation frequency is [7] 




(4) 
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Figure 1. shows the phase boundaries in the removal rate fi vs. diffusivity D space. On 
the left of the broken curve and above the continuous horizontal line the Turing state is 
stable. In what follows, we will concentrate on the region below the continuous horizontal 
hne in Fig.l where a Hopf state is available. 

Here we are going to consider the problem which is often treated in hydrodynamic 
instabilities - the one where the control parameter is given a sinusoidal temporal variation 
[8, 9, 10, 11]. To do so wc first linearize Eq.(l) around the fixed point B = = [1 + a)"^ 
and near the boundary /i = Thus we have 



(5) 



where the operator Lq reads 



d l-a 



L = { dt 1+a (l+(r)2 



(6) 



We note that the state obtained for /j, < is spatially homogeneous but temporally 
oscillating. We now endow jji with a small amplitude time modulation 

= /^o (1 + eCosiujt)) 

(7) 

and ask the question: what is the critical value of /Uq for which the spatially homogeneous 
oscillatory state is formed. 

The linearized equation now reads 

(8) 

The critical value of ^jlq can be expanded as 

A^o = A^oo + e/^oi + e^A*02 + •••• 

(9) 

where /ioo = and in Lq it is implied that /x = ^oo- In ^(1)) eigenvector for the 
homogeneous oscillatory state below the afore said Hopf-bifurcation boundary is 



^"^'^ - ( 2,..„a+.) )e-°* + c.c. 

\ /ioo+iti'o / 

(10) 



Thus we see that SAq and 5Bq has a constant phase difference since the production of B 
depends on A. At this point we consider this Phase difference 0(^) has an additive part 
which varies on a slower time scale er. So the structure of (f){fj,) is taken as 

(i) = (pc + S(P{t) 

(11) 
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where is the critical value of and can be easily obtained from Eq.(lO). The 8<p is 
expanded in powers of e as 

5(f) = 5(f)Q + e5(pi + e^502 + 

(12) 

Let us expand 5 A and 5B as 

5A — 5ao + eSai + e^Sa^ + .... 
5B = 5bo + e56i + e^5h2 + .... 

(13) 

and introduce the multiple time scales as 

t = tQ + er 

(14) 

where 

d d d 
dt dto dr 

(15) 

We now specialize to the case of parametric resonance i.e. to — 2u!q. The 0(e) equation is 

• 6ai \ f 

.-§^{5(f)o) + fj,oi[2{l + a)Sa+ - , ^ 

(16) 



Sbi J \6b+-§-{5(f)o) + //oi[2(l + a)Sa+ - b+] + ^[2(1 + a)5a_ - db^] 



On the right hand side of above equation, only the secular terms have been considered. 
The ± sign indicates ±iu! in the appropriate expressions of amplitudes Sao and 6bo of the 
0(0) solution. Now, the condition of occurrence of 0(e) solution i.e. the vanishing of 
secular terms lead to 



d .^^ , /xoi[2(l + a)5a+ - b+] + ^[2(1 + a)5a. - Sb. 
d-r^^"^'^ = Tb 



(17) 



Coming back to the original time scale and simplifying 

di^^"^'^ - - 2{^oo + 1) ) ~ 74^0 

(18) 

The condition that there will occur 0(e) solution is the occurrence of a mismatch in 
frequencies of oscillation of activator and inhibitor by the amount ^ — ^° ^ . 
One has to pull the Hopf boundary up by the amount 

fJ'OoifJ'OO — 1) 



2(/.oo + 1) 



(19) 
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to do away with such frequency mismatch between activator and inhibitor or in other 
words to get rid of weak disturbances. 

The solution of 0(e) equation for (a; = 2ijJq) will now definitely include e^'^^° term. 
This term will make secular terms appear in the next higher order. The 0(e) solution in 
g±i3a;o come as 

(20) 

The AstjQ being the determinant of Lq at frequency Scuq- Let us look at the O(e^) equation 
with only the secular part (i.e. the part of frequency ujq) on the right hand side. 



(5a2\^f 



V 5h2 J \6h+^mi) + //02[2(1 + a)5a+ - 6+] + ^[2(1 + a)5a^ - 5W] 



(21) 



Again the removal of secular terms results in 



d(l)i _ u>o 
dt /loo 



M A^02 + (— ) -T — - (— ) 



2 ^3u)o J 2 A3(jQ 



(22) 



So the condition for existence of O(e^) solution is a further mismatch of frequency of 
oscillation of the reacting species. A fine tuning of by 

//^00x2 1 
/^02 = (^) 



2 ' nlo + 9^i-f^ 

(23) 

can help getting rid of this frequency mismatch at this order too. In the above expression 
we have put 

(24) 

Thus we see that persistence to retain the frequency luq for the O(e^) solution causes 
a further shift in the Hopf-bifurcation boundary. 

The O(e^) solution will also have a e^*^'^°* part resulting from fii[2{l+a)Sdi — 6bi\ which 
will cause secular term to appear in the immediate higher order. The above mentioned 
part of the solution in O(e^) is of the form 

( S) - 2£kP(^ + -)^- - ^^"^] (.3.^.00 ) ^ 

(25) 
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This can be written in a simplified form as 



602 

5b2 



2A3a;o \^bi 



(26) 



Again the O(e^) solution will contain a e^iScuot part coming from the term //i[2(l+cr)5a2 — 
562] resulting in 
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and so on. In this way at 0(e)"' 

Sbr, 



2A3a;o 



{iSuJo + 1) 



5bi 



(27) 



2A3a;o 



{iSuJo + 1) 



(n-l) 



6ai 
Sbi 



(28) 



Thus we get a nontrivial flow of secular term generators at all orders producing phase 
instabilities. Such a band of phase instability is generated as a result of sub-harmonic 
response to the external force. These instabilities would possibly show up in weak phase 
turbulence under forcing. 

In the conclusion we would like to mention that, a multiple scale perturbation anal- 
ysis of forced Hopf state reveals the presence of instabihties which are responsible for a 
slow frequency drift between the two reacting species. We would like to focus on the 
point of instantaneous varied response of activator and inhibitors is the basic cause of 
generation of such instabilities. This type of differential response to an applied force can 
always occur when one of the reactants depends on the other for its production and thus 
allowing for a delay. A persistence of the reacting species for oscillating in unison can 
cause in a slow overall phase drift of the system. Such oscillatory regions separated by 
continuous distribution of relative phases has been experimentally observed by Lin et al. 
in the forced Belousov-Zhabotinsky system at low forcing amplitude [4]. At high enough 
forcing amplitude they have got well defined n phase separated 2:1 resonant pattern. It is 
important to note that a gradual increase in the parameter in our analysis, to get rid of 
spurious instabilities is also the same as enhancing the forcing amplitude. Consideration 
of large scale spatial phase variations under resonant forcing can easily be shown to result 
in an inhomogeneous diffusion equation in with a complex inhomogeneity C originating 
as a result of external forcing. This type of a situation can result in spatial instabilities in 
one of the reacting species. Thus other spatio-temporal instabilities at larger scales can 
also be shown to exist, which influences the situation under forcing. 
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Figure caption 

Figure 1 shows the phase diagram of the one dimentional Gierer-Meinhardt model on 
removal rate fj, vs. diffusivity D space. The continuous line is a Hopf bifurcation bound- 
ary whereas the broken hne separates steady Turing state (on the left) from the basic 
homogeneous steady state. 
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